Dynamic density functional study of a driven colloidal particle in polymer solutions 
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The Dynamic Density Functional (DDF) theory and standard Brownian dynamics simulations 
(BDS) are used to study the drifting effects of a coUoidal particle in a polymer solution, both for 
ideal and interacting polymers. The structure of the stationary density distributions and the total 
induced current are analyzed for different drifting rates. We find good agreement with the BDS, 
which gives support to the assumptions of the DDF theory. The qualitative aspect of the density 
distribution are discussed and compared to recent results for driven colloids in one-dimensional 
channels and to analytical expansions for the ideal solution limit. 
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INTRODUCTION 

Mixtures of colloids and non- adsorbing polymer coils 
have attracted much attention over the past decade both 
experimentally and theoretically fl|. They provide an 
excellent model system to understand the generic equi- 
librium and nonequilibrium physics of multicomponent 
colloidal mixtures due to the fact that the interactions 
between the constituents can be tailored The theo- 
retical approaches to the study of the equilibrium phase 
behavior of colloid-polymer mixtures were largely based 
on the Asakura-Oosawa (AO) model, in which the 
chains are modeled as ideal particles experiencing a hard- 
core repulsion with the colloids. Recently, extensive com- 
puter simulations Q revealed qualitative differences in 
the phase diagrams when interactions between the poly- 
mers are included. However, the off-equilibrium behavior 
of these systems is far from being understood. In this 
paper, we present results based on a recently proposed 
dynamical density functional (DDF) formalism [5|,|6j and 
we demonstrate that the latter is capable of describing 
out-of-equilibrium diffusive processes at the Brownian 
timescale. The advantage of the DDF theory is the fact 
that particle interactions are included once a good ap- 
proximation for the equilibrium functional is known and 
it is well suited to treat different external potentials. 

The system under consideration is a colloidal parti- 
cle being dragged at a constant rate c (e.g. by gravita- 
tion, electric or magnetic fields, or by optical clamps 0) 
through a solution of polymers in a light solvent. The 
spherical colloid is represented by an external potential 
V^xt(r,t) = T4xt(|r'|), where r' = r — ctz, is the coordi- 
nates in the reference framework of the colloidal particle. 
The solvent provides the rest-framework for the Langevin 
dynamics of the polymers, which have a mobility Fq, 
connected through the Einstein relation to the diffusion 
constant and the inverse thermal energy, /3 — [ksT)^^ . 
Assuming that the polymer gyration radius and the col- 



loidal particle have similar size, a ~ lO^^m, and with the 
viscosity of a typical solvent at room temperature, the 
natural units for the shifting rate Fq/ (/3cr) would be in 
the range of 10~*m/s, and we may safely neglect the hy- 
drodynamic effects of the light solvent. The much heavier 
globular polymers would feel the competing effects of the 
solvent, as rest framework for their Brownian dynamics, 
and the moving colloidal particle drifts with respect to 
the polymers at a constant rate (drift velocity) c. The 
deterministic Dynamic Density Functional (DDF) theory 
0,0, is an extension of the Density Functional (DF) for- 
malism to off-equilibrium systems, which includes exactly 
the ideal gas and the external potential contributions to 
the free energy, and it represents the correlations out 
of equilibrium by those of an equilibrium system with 
the same density distribution. With this hypothesis, and 
the interpretation of the density p{r,t) as the average 
of the instantaneous density over the random noise in 
the molecular Langevin dynamics, the theory enables the 
use of the well developed approximations for the equilib- 
rium Hclmholtz free energy in DF theory, which are usu- 
ally split into the ideal gas and interaction contributions, 
J^[p] = ^ideai[p] + AJ^Ip]. The central DDF equation for 
the time-dependent density distribution is 
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From any initial distribution of polymers, the time evo- 
lution would take the system towards to a stationary den- 
sity distribution p{r,t) — p{r — zct) = p(r'), shifting at 
the same rate c as the external potential. This distribu- 
tion is the most relevant property of the system and it 
corresponds to the solution of the functional equation 



= 0, (2) 



where the time dependence is fully adsorbed into the co- 
ordinate r'. The solution of Eq.Q has been explored 
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for systems with one-dimensional (ID) dependence of the 
potential barrier Kxt(r') = Voxtiz'), along the shift- 
ing direction; the similarities and the differences with the 
Euler-Lagrange equation for the DF theory of equilibrium 
systems were analyzed there both for ideal and interact- 
ing systems. However, the three-dimensional (3D) geom- 
etry of the present problem requires a different analysis, 
since the zero-divergence requirement for the brackets in 
Eq.(|2Il leaves open a much wider functional space in 3D 
than in ID. 

NON-INTERACTING POLYMERS 

Starting in the spirit of the AO model, for non- 
interacting polymers, AT = 0, Eq. (0) becomes a linear 
Fokker-Planck equation 

VV + V(p- V/3\4) = 0, (3) 

with a "kinetic potential" given by /3Vfe(r') = /3T4xt(r') + 
cz', as a function of the reduced shifting rate c = Pc/Tq 
with inverse length units. In Fig. l(a,b) we present nu- 
merical solutions for the density distribution of the ideal 
case with bulk density poa^ = 1, under the effects of the 
external potential 

K=xt(r') = Vo exp(-|r7a|6), (4) 

with /3Vq — 10, to represent the soft repulsion between 
the polymers and the shifting colloidal particle. The bulk 
density po is the value of the density far away from the 
external potential and is obviously the same both for 
the equilibrium and the nonequilibrium, driven system; 
moreover, for the ideal non-interacting system po just 
provides an arbitrary factor to /o(r')- 

The density distribution around the external potential 
/o(r') has axial symmetry and exhibits a cap- like structure 
with p(r') > Po in the front {z' > cr), formed by the 
polymers being pushed by the moving repulsive external 
potential. These particles escape around the colloidal 
particle creating a skirt for z' < —a, which, together 
with the hole ( p(r') < po) left behind by the potential, 
form a wake structure extending much further than the 
front structure for z' ^ a. An increased shifting rate 
( cfj = 10 in Fig. lb compared to ccr = 1 in Fig. la), 
enhances these characteristics, with a higher density front 
and a narrower skirt reaching further away behind the 
shifting colloidal particle. 

The transverse integral of /o(r') — po over x'y' plane 
as function of z' gives an excess of polymers in the 
front region but it vanishes in the wake region, as the 
positive wings exactly compensate the depletion close 
to the z'-axis. This is reminiscent of the ID result 
for a shifting repulsive barrier with a front density 
p{z') = Po + Ae'xp{—cz') and no wake , p{z') = po, be- 
hind. The absence of excess molecules in the wake seems 



to be a generic characteristics of stationary states under 
constant bulk boundary conditions with purely relaxative 
dynamics of the ideal gas molecules. Using cylindrical co- 
ordinates r' = (i?, 0', z') we now explore analytically the 
asymptotic forms of both the front and the wake regions, 
where Kxt(?'') = reduces Eq.® to: 

1 d (^dp{v')\ , d^pjv') , _dp{v') 

RdR r ^ j + + = ^- 

Through a Hankel transform, the solution for the cylin- 
drically symmetric p(r') is, 

p{R,z')^po+ daaf{a)Jo[aR]e-'^''] (6) 

^0 



where /3± = | ± y/(|)2 + 0;^ and Jq is the zeroth or- 
der Bessel function. Far away from the external poten- 
tial, a fast convergence of the Hankel components f{a) 

is observed, and the relevant features come from their 
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behavior for a <^c. Hence we may use /?+ « c -f ^ and 

2 

/3_ « — for the front and the wake regions respec- 
tively. The expansion of f±{a) as an even polynomial 
function for small a and the zero wake requirement lead 
to f-{a) Aio^ + A2a^ + ... in the wake , while at the 
front we expect /+(«) « -Bq + B^a^ + B2a^ + ... Thus 

p{R,z')^Po + e-^^- e--[Bo^ 
+B,^ {a' ~ 4) + (a^ - I6a' + 32) + ...], (7) 

for z' 3> cr, and 

„2 y;4 

p{R,z')^po + e--[Ai—{a^-A) 

+^2^(a'-16a2 + 32) + ...]; (8) 

for z' <C — cr; where w — yjc/ \ z' \ and a — wR. 
Although the amplitudes of these contributions depend 
on the particular external potential, the asymptotic de- 
cay forms are generic. For a fixed z', the structure in the 
transverse plane is given by a gaussian times a polyno- 
mial function. In our case the positive values of Ai and 
Bo lead to a maximum front density at i? = 0, while the 
leading term at the wake gives a minimum at i? = and 
maximum at i? = 2^/^y^| z' \ /c, producing the cap-like 
structure with parabolic shape shown in Fig. l(a,b). Be- 
sides the lack of Ao, the qualitative difference between 
the advancing front and the wake comes from the ex- 
ponential decay exp(— cz') behavior of the front , which 
restricts the density structure to the neighborhood of the 
external potential, as in the ID result, while for fixed a 
the wake structure decays as inverse powers of z', with 
a 1/z'^ leading term. We have tested these analytical 
predictions with the numerical solution of Eq.(|2Jl, and 
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got quantitative agreement for z' > 1.25cr (ccr = 10) and 
z' > 3(7 (ccr = 1) in the front region, and z' < —7a in 
the wake region for any c. The contributions from higher 
order terms in Eq.® appear to be more important than 
in Eq.Q, although the quahtative aspect of the wake is 
already well represented by the first term in Eq.®. 

Nevertheless, we have to point that, contrary to the ID 
case 1^ , the front and the wake regions in our 3D system 
are in fact connected through the regions, with small \z'\ 
but large R, were the external potential created by the 
shifting colloidal particle on the polymers vanishes; so 
that the solution of Eq.Q should have a unique analytic 
form, common to the front and the wake region. This 
has the obvious difficulty of using both the positive (/3+) 
and the negative (/3_) decay constants for the positive 
(front) and negative (wake) values of z' , leading to ex- 
ponential growth of their respective contributions, which 
may only be canceled by the appropriate behaviour of 
/(a) for large a, beyond the Taylor expansion used in 
Ea. (|7l8(l . The good comparison of our numerical solu- 
tions with the analytic results Eq.(|Zj) for the front and 
Eq.® for the wake reflects a local asymptotic conver- 
gence which is quite useful to understand the qualitative 
features of the stationary density distribution, but which 
cannot be taken as an exact global asymptotic result. It 
has to be pointed that the use of spherical, rather than 
cylindrical coordinates to solve Eq. © also leads to prob- 
lems of convergence, as the parabolic structure of the 
wake implies the entanglement of the radial and the an- 
gular coordinates. 
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FIG. 1: Steady state contour density field of ideal (a,b) and 
interacting (c,d) polymers created by a driven colloidal par- 
ticle. Shown is the density in the x'z'-plane where the center 
of the colloid is located. The colloid is moved along the z'- 
axis at a velocity ac = 1 (a,c) and ac — 10 (b,d). Bright 
regions correspond to low densities while dark regions show 
high densities. 



THE EFFECT OF POLYMER INTERACTION 

For the case of interacting polymers the steric effects 
lead to an effective repulsion between them, which we 
model by the ultra-soft gaussian pair potential 0, 0, 

ElEllillIi 

0(ry) = C/oexp(-4/a2), (9) 

where rij is the interparticle distance and (3Uo = 1. Both 
for equilibrium [l3 . Il^ ITsf and dynamical properties 
^3 the excess free energy density functional of this model 
has been successfully approximated by a pure mean field 
(or random- phase approximation (RPA)) form 

= ^ /^'^/ d'r'Hlr - r'l) p(r)p(r'). (10) 

Instead of solving the integro-differential Eq.® with 
this model, we have obtained the steady state distribu- 
tion by the time integration of Eq.(^ from a uniform 
density initial state. The stationary structure around 
the colloidal particle are reached with short integration 
times, as presented in Fig. 1 for po<J^ = 1 (which rep- 
resents a fairly dense polymer solution) and c = Icr (c) 



and c = lOcr (d). At the low velocity ccr = 1 the influ- 
ence of interactions between the polymers is very strong. 
The spherical layering structure created around the col- 
loidal particle by the polymer-polymer repulsion is much 
stronger than the front - wake asymmetry induced by 
the dragged colloid; the extension of the wake behind the 
moving particle is strongly reduced by the much lower 
bulk osmotic compressibility of the interacting system, 
which facilitates the filling of the axial hole by radial 
currents. At the higher ccr — 10 shifting rate in Fig. 1(d) 
the effects of the interactions are much weaker. Although 
there is still a clear shortening of the wake , explained by 
the lower osmotic compressibility, the main qualitative 
change with respect to the ideal solution result in Fig. 
1(b) is that the layering created by the packing effects 
produces a double cup-like structure, reaching further 
away from the z'-axis. 

In Fig. 2 we present a quantitative view of our results; 
we plot the polymer density as a function of the distance 
R to the z'-axis for fixed values of z'/cr = ±1. The solu- 
tion of the DDF approach is compared to standard Brow- 
nian dynamics simulation (BDS) 0|. Here, the stochas- 
tic Langevin equations for the overdamped colloidal mo- 
tion of N particles with trajectories ri{t) {i = 1,...,N) 



read as 

ro'^ = -Vr. 5^0(1 r.-r, |)+Fe.tW+ro"'c+Ff)(t). 

(11) 

There are different forces acting onto the cohoidal par- 
ticles: first there is the force attributed to inter-particle 
interactions, secondly there is the external field Foxt due 
to the colloidal particle, r,^^c is the driving force, and 

fry \ 

finally the random forces describe the kicks of the 
solvent molecules acting onto the ith colloidal particle. 
These kicks are Gaussian random numbers with zero 
mean, F^^^* = 0, and variance 



(FrOa(i)(Fr)/5(i') = 2kBTToSa.f3S,,S{t - t') (12) 

The subscripts a and /3 stand for the three Cartesian 
components. The simulations were carried out with 
N — 1000 particles, periodic boundary conditions in all 
directions and box sizes = Ly = 8a and = 16cr. 
After an equilibration time of 10^ timesteps, statistics 
were gathered over a period of 2 ■ 10^ timesteps. In both, 
simulation and numerical solution of the DDF, the den- 
sity is averaged over rings at z' and radius R with a cross 
section cr^/4. 

The good agreement between the BDS data and the 
DDF results gives support to both the mean field approxi- 
mation H1U|I used to describe the effect of the interactions, 
and to the DDF borrowing of the equilibrium intermolec- 
ular forces, as functionals of the instantaneous density 
distribution |^ Q. The comparisons between the ideal 
solution and the interacting system shows again that the 
polymer interactions affect much more the features of the 
density distribution for low c than for high c. For ca = 1 
the polymer layer around the colloidal shows only a slight 
asymmetry for cr ^ z' ^ —a, while the kinetic effects in 
ideal gas create a maximum at the front z' ^ a. 

Qualitatively we may associate the behavior of the 
front structure to the direct kinetic effect of the advancing 
spherical repulsive potential created by the colloidal par- 
ticle. As c grows the kinetic constrain on their brownian 
trajectories becomes the dominant factor for the move- 
ment of the polymers at the front . The polymer-polymer 
interaction plays then a minor role, so that for ca ^ 1 
the main peak in the front structure becomes similar for 
ideal and interacting polymers. The structure of the wake 
is determined by the diffusion from the bulk solution to 
fill the void left behind the colloidal particle, the effect of 
the higher osmotic pressure accelerates that process and 
produces a weaker wake than in the ideal solution limit. 
On the opposite extreme, for very low shifting rates of 
the colloidal particle, the effects of the polymer-polymer 
interactions are very important, both at the front and at 
the wake structures. At high polymer concentrations the 
structure around the colloidal particle is dominated by 
the steric repulsions between polymers. The relatively 
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FIG. 2: Steady state density profiles of the polymers around 
a driven colloid at z'/ct = plotted as a function of the radial 
distance R/a from the z'-axis for fixed z' /a = ±1. The shift- 
ing rates are (a) ac = 1 and (b) crc = 10. Dashed lines are 
the results for ideal polymers, while the curves for interacting 
polymers are plotted with solid lines. The symbols (circles) 
are the simulation results. 

rigid structure of molecular layers is shifted, with little 
deformations, by the moving colloidal particle. 

DISCUSSION 

As the first point in our discussion, we consider the 
total excess of polymers AA^, produced by K;xt(r') over 
the uniform bulk density. This is a relevant data since 
cAA^ is the total polymer current, which requires a total 
force cFqAA^ provided by the colloid on the polymers. 
There is a generic DDF relation 



AN = 




for the stationary density distributions, which allows 
to calculate the total excess from the density distribu- 
tion in the neighborhood of the external potential. The 
results for the systems in Fig(l) are: AA"a = 3.05, 
AA^f, = 1.883, AAc = 2.53 and ATV^ = 1.879. Again, 
the difference between the interacting and the ideal cases 
is strongly reduced at ca = 10, with respect to ccr = 1. 
It is remarkable that the static equilibrium result for 
AA^ in the ideal gas case would be negative (as p(r) = 
Po exp(— /3T4xt(r)) < Po), but the stationary excess AA 
is positive and grows as c decreases. This may be un- 
derstood from the analytical ID result ^ for the front 
p{z') — po ~ c exp(— cz'), which vanishes locally as c — > 0, 
but it still gives a positive integral over-compensates 
the depletion inside the potential barrier and produces 
AA^ > 0, consistently with the sign (and value) of the to- 
tal force. The difference between the equilibrium (c = 0) 
density distribution and that of a stationary state at ar- 
bitrarily small but positive c is remarkable. The appar- 
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ent paradox come from the concept of stationary state, 
which would appear after a short transient period when 
c is large, but it would require diverging times as c — + 0. 
The very weak but extended structure of the exponen- 
tial front in p{r') for ca <^ 1 would never be observed in 
practice, and the transient states p{r, t) for any reason- 
able t would be very similar to the equilibrium structure 
for c = 0. Nevertheless, the strongly anisotropic den- 
sity distributions, with non-trivial global effects even for 
very low c, suggest important effects on the interaction 
between two driven colloids in a bath of quiescent Brow- 
nian particles qualitatively different from effective 
interactions in equilibrium [l|, [iflj . 

Finally we comment on the relevance of the boundary 
conditions and the system dimension by comparing the 
present results with the ID system explored in a pre- 
vious work 1^. Obviously, in the system explored here 
the effects on the bulk polymer solution are limited to 
the neighbourhood of the single colloidal particle. If we 
consider a finite concentration of colloidal particles, all 
being drifted at the same rate c with respect to the sta- 
tionary framework of the solvent, there would be a finite 
induced polymer current per unit volume. The 3D struc- 
ture, which offers easy paths for the polymers to escape 
from the colloidal particles, would probably make unfea- 
sible the approach to the full-drift regime discussed for 
ID systems, in which nearly all the particles move along 
the shifting potential. Another possible problem in 3D 
which may be of interest are currents through a struc- 
tured barrier with holes or slits. However, we have to be 
aware of the intrinsic limitations of our DDF approach, 
particularly in the treatment of the solvent as an inert ref- 
erence framework for the Langevin dynamics of the poly- 
mers, which is not affected by the shifting external poten- 
tial (or colloidal particle). Altogether, we may conclude 
that the DDF offers a good theoretical tool to explore dy- 
namical problems in polymers solutions subject to time 
dependent external potentials which do not affect the sol- 
vent, and the formalism may also be used to predict, in 
good agreement with BDS, the density structures cre- 
ated around colloidal particles or similar sized molecules 
moving slowly with respect to the solvent. However, the 
extension to problems in which the solvent plays a more 
direct role should be regarded with caution, as they may 



require a more symmetrical treatment of the solute and 
the solvent. 
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